Dissipation-driven quantum phase transition in superconductor-graphene systems 
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We show that a system of Josephson junctions coupled via low-resistance tunneling contacts to graphene 
substrate(s) may effectively operate as a current switching device. The effect is based on the dissipation-driven 
superconductor-to-insulator quantum phase transition, which happens due to the interplay of the Josephson 
effect and Coulomb blockade. Coupling to a graphene substrate with gapless excitations further enhances charge 
fluctuations favoring superconductivity. The effect is shown to scale exponentially with the Fermi energy in 
graphene, which can be controlled by the gate voltage. We develop a theory, which quantitatively describes 
the quantum phase transition in a two-dimensional Josephson junction array, but it is expected to provide a 
reliable qualitative description for one-dimensional systems as well. We argue that the local effect of dissipation- 
induced enhancement of superconductivity is very robust and a similar sharp crossover should be present in finite 
Josephson junction chains. 



Artificially fabricated Josephson-junction arrays (JJA) offer 
a unique playground for studying quantum phase transitions 
(QPT) [1]. The transitions in JJA occur due to the competi- 
tion between the Josephson coupling, E j, which favors a glob- 
ally ordered state, and the charging energy, E c , which leads 
to Coulomb blockade of Cooper-pair tunneling and enhances 
quantum fluctuations of the superconducting (SC) phase. At 
zero temperature, the QPT from a globally superconducting 
to an insulating phase EL fl IS El 0] occurs, roughly, when 
the Josephson energy Ej becomes smaller than the charging 
energy. Another key factor in determining the ground state is 
dissipation, which is present whenever the SC system is con- 
nected to a reservoir of gapless single-particle excitations 10]. 
The main effect of the dissipation in JJA is a suppression of 
quantum phase fluctuations. Taking into account the phase- 
charge uncertainty relation, the dissipation enhances fluctu- 
ations of the charge and, hence, stabilizes the SC phase||8|]. 
This type of dissipative QPT has been considered previously 
by Feigelman and Larkin [9] in the model of a regular 2D 
proximity-coupled JJA and by Galitski and Larkin fl 1 011 in a 
disorder-induced random Josephson network. In both cases, it 
was found that the effect of dissipation on the transition point 
is exponential, i.e., the critical Josephson coupling scales ex- 
ponentially with the Andreev conductance. The Andreev con- 
ductance and hence the degree to which the stabilizing effect 
of dissipation is important obviously depend on the density 
of states of gapless excitations. Thus, by controlling the lat- 
ter, one can tune transitions between a global superconductor 
and an insulator. This observation provides strong motivation 
for studying superconductors coupled to a graphene substrate, 
where the density of carriers can be easily tuned by gate volt- 
age from essentially zero at the Dirac point (no "Ohmic" dis- 
sipation) to very large values in the electron-doped graphene 
(strong dissipation). 

In this Letter, we propose to study the superconductor- 
insulator phase transition (SIT) in a Josephson-junction ar- 
ray in a tunneling contact with a graphene layer (or layers) 
(see Fig. [U, which acts as a source of gapless quasiparticles. 



The graphene substrate provides a unique possibility to con- 
trol the dissipation strength via the gate voltage and thereby 
tune the dissipation-driven QPT. Hence, the system may be 
used as a current switching device. While the physics of the 
underlying effect is intuitively quite clear, the formal descrip- 
tion of the transition developed in this Letter is technically 
non-trivial: First, we use the tunneling Hamiltonian formal- 
ism and elements of random matrix theory to derive the ef- 
fective phase fluctuation action of a small SC grain coupled 
to a graphene substrate. The dissipation kernel, K(r), shows 
a crossover from the Ohmic dissipation behavior Kir) oc r~ 2 
in the electron-doped graphene to extremely weak dissipation 
K(t) oc t -4 at the Dirac point. Second, we develop a mean- 
field theory of the SIT and show that the quantum critical point 
is determined by the single-grain phase correlator, which is to 
be calculated using the dissipative effective action. To calcu- 
late the phase correlator we use the two-loop renormalization 
group (RG) results from a related spin model and determine a 
critical voltage, V C9 at which the transition occurs: For V > V c , 
the system is a superconductor; for V < V c , the system is an 
insulator. 

Our theoretical model is an array of SC grains connected 
with each other with the Josephson junctions and connected 
via tunnel contacts to a graphene substrate, see Fig. [T] For 
T — > 0, one can neglect massive fluctuations of the amplitude 
of the order parameter A in the grain, and describe the dynam- 
ics of the system in terms of the phase-only imaginary-time 
effective action {% = 1) 



(1) 



where (p t is the phase of the order parameter on the i-th grain. 
Here, for simplicity, we assume that the Josephson and charg- 
ing energy are the same for all grains. However, this assump- 
tion is not essential for our results. 

We now consider the effect of graphene gapless excitations 
on the phase coherence of a single SC grain. We study here 
the situation when the SC grain of radius R lies on top of the 



2 




^ CO 



FIG. 1: (Color online) a) A Josephson-junction array on top of a 
graphene sheet is shown. The grains are coupled by the Josephson 
junctions with the coupling strength Ej. The distance between the 
grains L is assumed to be large compared to the SC coherence length 
Thus, the coherent transport through the graphene is neglected 
here. Alternatively, the SC grains in the array can be coupled to dif- 
ferent graphene sheets. Then, there is no coherent transport through 
the graphene due to physical separation between the sheets, b) The 
superconductor/graphene interface. The graphene and superconduc- 
tor are connected through the tunnel barrier. The chemical potential 
in graphene can be tuned with the gate voltage V g . 



graphene sheet, see Fig. \T\ In this planar geometry the tun- 
nel junction does not break internal symmetries of graphene, 
and thus does not modify the spectrum of the excitations. In 
the limit of low transparency tunnel barrier, the transport be- 
tween superconductor and graphene can be described by the 
tunneling Hamiltonian 



H T = t^fd z Y J f A d 2 r [^(r)^V)+H.c.] , 



(2) 



where *F^(r) and ^(r) are the electron operators in 
graphene and superconductor, respectively. Here d z is the 
thickness of the grain, and t and A are the tunneling matrix 
element and the area of the junction, respectively. 

Assuming that the SC gap energy is sufficiently large, the 
main contribution to the subgap transport originates from An- 
dreev processes, which involve correlated tunneling of two 
electrons from/to the graphene. In the fourth order perturba- 
tion theory in tunneling t (see Fig. |2]), the contribution of the 
Andreev processes to the dynamics of the phase is given by 
the following effective action 



S A = 24fd: 



Ja i=l.A 



(3) 



Re [f*(x u x 2 )F(x 3 , x 4 )G (g \x u x 3 )G (g \x 2 , x 4 )] , 

with x = {r, r}. Here G^ g \x, x') and F(x; x') are imaginary- 
time Green's functions for graphene and superconductor, re- 
spectively. The latter is defined as 

F*(ri , r 2 ; ti , r 2 ) = -e«nMT 2 )]/2 J^ Xn ( ri ) Xn (r 2 )u n v n 

n 

[0(n -t 2 )^- (Ti - T2 ^0(^)-0(t 2 -t 1 )^- (Ti - T2) ^0(-^)] . 



Here E n = ^js^ + A 2 , u n and v n are Bogoliubov coherence 
factors u^/vl = (1 + s„/E n )/2, and s n and;^(r) are the eigen- 
values and eigenfunctions of the single particle Hamiltonian 
of the grain, which includes random potential due to impuri- 
ties and boundaries of the island. 

In order to derive low-energy effective action S a due to An- 
dreev processes, it is convenient to separate the fast and slow 
times k and r: t\ = r + k/2 and r 2 = r - k/2. Since the su- 
perconducting Green's function decays exponentially on the 
time scale of order A -1 , one can neglect k in the dynamics of 
the phase, which evolves on much longer time scales. Tak- 
ing advantage of this separation of scales and performing the 
integral over the fast time scales, we obtain the following ex- 
pression for Sa 



S A ~24t 4 d 2 



\ [I dridrdr' cos[(f>(T)-(f)(T')] 

Ja iii4 

xF u F 34 G^(t-t')G^(t'-t). 



(4) 



Here F t j is the anomalous Green's function at zero frequency, 
i.e Fij = - XnXn(rdXn(rj)A/El. 

It is well-known that Andreev transport is sensitive to disor- 
der fllll . Therefore, in order to calculate the effective action, 
one has to take into account spatial correlations 
and average the action over the random realization of the 
wavefunctions in the SC grain resulting from the scattering of 
electrons from the grain boundaries and impurities. We per- 
form this averaging using exact eigenstates technique assum- 
ing that the grain is sufficiently small. Our approach accounts 
for the enhancement of the tunneling rate due to the coherent 
back- scattering of electrons to the tunnel junction. The corre- 
lation function (F12F34) in the grain consists of reducible and 
irreducible parts: 

<^12^34> * <Fi2><F 3 4> + (Fl2F 34 ) ir . (5) 

The reducible part can be easily calculated 

/Z7 \ 71 t ^ sinfeln - r 2 |) Jirpl ^ 

(Fn) = -v F /i2 with f u = ——. —e * , (6) 

2 k F \ri - r 2 \ 

where hp, I and vp are the Fermi wave vector, mean free path, 
and density of states at the Fermi level in the grain, respec- 
tively. 

The irreducible correlation function is calculated assuming 
that SC grains are small, and the Thouless energy E T is the 
largest relevant energy scale in the problem E T ~ D/R 2 » 
A, E c ,Ej, where D is the diffusion constant, R is the radius of 
the grain, and 6 is the level spacing in the grain. In this limit, 
the electron diffusion time in the grain r D = 1 /E T is much 
smaller than the time the system dwells in the virtual state 
with one unpaired electron ~ 1 / A. Since an electron in the 
virtual state covers the entire available phase space, one can 
calculate the irreducible correlation function within Random 
Matrix Theory [14] finding that it acquires the universal form 

(FnF 34 )ir = ~ 4 ^(/i4/23 + /13/24). Combining this equation 
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FIG. 2: The diagrams describing the correlated two-electron tunnel- 
ing process (Andreev process) from/to the SC grain are shown. 



with Eq. ©, we obtain 

(^12^34 > = ^V 2 F fnf34 + ~V* F (fl4fa + fafa). CO 

Using Eqs. (|7]) and ©, and carrying out the spatial integrals, 
we find 



6ti f 
SAkp 



J drdr'c 



COS[0(T)-0(T / )]G<f(T-T / )G^(T / -T). (8) 



This term originates from the integration of the irreducible 
correlation function and is parametrically larger than the re- 
ducible term by the factor of k^AS/A » 1 lllllll3ll . After sub- 
stituting the expressions for the graphene Green's functions 
(see, e.g., Ref. 1 15]) and assuming that [i > 0, we find 

S A *-G J iTdT'K(T-T')C0s[(P(T)-<P(T')] , G = - , (9) 

where y is the Fermi- velocity in graphene. The kernel K(r) is 
described by the following function 

K(t) = ^^{0ur) 2 + 2[l- y u|r|]^ |T| - l), (10) 

which exhibits a crossover from r~ 4 at the Dirac point to the 
usual Ohmic behavior r~ 2 at fi\r\ ~ 1. For realistic exper- 
imental parameters the chemical potential [i is much larger 
than the characteristic timescale for the dynamics of the phase 
r < E~ l , i.e. jj » E c . Thus, the asymptotic form of the 
effective action becomes 



, cos [0(t)-0(tO] 
drdr : — , rj 



(t-tO 



./\2 



k 2 F 7 2 



(11) 



The important difference between action (ITTI) and the dissipa- 
tive action describing resistively- shunted Josephson junction, 
is that action (fTTb is "compact" or periodic in the phase dif- 
ference, and thus correctly describes the fact that the charge 
on the grain is quantized and can change by 2e only. From 
Eq. (fTTb , we see that graphene as a source of gapless ex- 
citations provides the possibility to change the dissipation 
strength directly by changing the chemical potential, which 
is tied to the gate voltage. 

Dissipation-driven QPT. Combining Eqs. (Q]) and (fTTb , we 
obtain the full action for the system: S =Sc + Sj + Sa- To 



derive effective action describing the transition in the JJA, we 
first write a partition function in path integral representation 
and then use Hubbard- Stratonovich transformation to decou- 
ple Josephson term by introducing an auxiliary field if/i(r) cou- 
pled linearly to e m<J \ Then, the partition function becomes 
Z = Zo f T)ifsexp(-S [if/]), where S [if/] is given by 



-ln^expU J dr^e^Vi 



(12) 



(r) + c.c. 



Here we introduced the symmetric matrix w^, which describes 
Josephson tunneling: matrix elements of Wij are equal to Ej 
for nearest neighbors and zero otherwise. The expectation 
value in Eq. (IT2l) is taken with respect to the single-site ac- 
tion So = Sc + Sa- To study QPT at the mean field level, 
one may perform cumulant expansion of the second term in 
the action S [if/] in powers of if/ and arrive at the effective com- 
plex iff 4 field theory L16J . The phase boundary between macro- 
scopically superconducting and insulating state of JJA can be 
obtained by setting the coefficient r in front of \if/\ 2 to zero: 



r oc 



^-If drie^-^h = 0. (13) 



Here z is the coordination number of the lattice and averag- 
ing is taken with respect to So. When calculating the cor- 
relation function we assume that the dissipation strength n 
is large. Then, the second term in So dominates at low fre- 
quencies and one can neglect the influence of the charging 
energy term Sc, which serves as the ultra-violet cut-off. Un- 
der these conditions, the correlation function can be mapped 
on the long-range ferromagnetic spin chain first 
considered by Kosterlitz JI17I1 . Later, the critical behavior 
was studied in Ref. 1191 [20|], where the asymptotic behav- 
ior of the spin- spin correlation function was obtained using 
two-loop renormalization group (RG). Adopting the results of 
Refs. [y,ll7lll9ll20|] to our problem, we get 



(e 



i(f(r)-i(p(0) 



>0 



(?); 



A 1 «T«T C 
T»T C , 



(14) 



where A is the ultraviolet cut-off, A ~ 2nE c r], and r c is the 
correlation time, which can be calculated using the RG for 
n. Since the compact dissipation term proportional to n is 
not Gaussian, it gets renormalized when integrating out high- 
frequency degrees of freedom resulting in the following flow 
equations 



drj 



dln(AT) 



1 

'2t? 



1 



(27T 2 ) 2 n 



(15) 



where the rhs is the beginning of a Taylor expansion in 1 /rj. 
By integrating this equation between the initial value 7/(0) = n 
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FIG. 3: (Color online). Phase diagram for the SIT is shown. By 
tuning the applied gate voltage V g , one can change the density of 
states of gapless excitations in graphene, and thus drive the phase 
transition between globally superconducting and insulating states of 
JJA. The solid (black) line represents the region of the validity of RG. 



and final value 27r 2 77(ln[Ar c ]) 
lation time as 



1 , we can estimate the corre- 



(4n 3 E c r] 2 )- 1 exp (2^77). 



(16) 




It is assumed here that 2n 2 r] » 1. Using Eqs. (ITU) and (IT6l) 
and carrying out the time integral, we find expression defining 
the phase boundary C z Ej = 47T 3 E c nl exp (-27t 2 t/ c ), where C z 
is a numerical constant depending on the coordination number 
z. The rhs of the above equation corresponds to the effective 
charging energy of the grain E* renormalized by the dissipa- 
tion [2l[ H2, 12311 . The critical dissipation strength, n c , and 
the corresponding critical voltage, V c , at which the QPT oc- 
curs can be expressed via Lambert-W function by inverting 
the transcendent algebraic equation above: 



(17) 



Here Vo = j^j= with hp, J and G being the Fermi-momentum 
in the SC grains, the Fermi-velocity in graphene, and the di- 
mensionless constant defining the transparency of the tunnel 
barrier (see Eq. ©), respectively. The phase boundary be- 
tween globally superconducting and insulating phases of JJA 
is shown in Fig. [3] For the theory based on the RG procedure 
to be formally valid, we need ln(E c /Ej) » 1. This, how- 
ever, is a mathematical rather than physical constraint and for 
all practical purposes (i.e., experiment) what is important is 
the existence of the transition itself; the applicability of RG 
methods and the exact location of the non-universal "critical 
voltage" are not essential. A sharp transition should certainly 
exist if E c » Ej and perhaps even for E c > Ej. In the oppo- 
site limit the system is already superconducting without any 
substrate and there is no dissipation-driven effect. Since the 
technology for making SC grains with required E c /E j is well- 
developed, the conditions for the observation of the QPT are 
certainly experimentally feasible. 



We also emphasize here that since the dissipation-driven 
transition is intrinsically local, it is expected to survive in 
one-dimensional chains (where the mean-field theory breaks 
down and the QPT is of Kosterlitz-Thouless type, see, e.g., 
Ref. [24]). Again, the exact location of the transition point 
and the critical behavior would be different, but the effect 
itself should be present. Moreover, the same local argu- 
ment suggests that a sharp voltage-induced crossover (cf. 
with the model of a shunt resistor coupled locally to the SC 
grain |25ll26|]) in the IV-curves should be present even in finite 
chains proximity coupled to graphene, similar to those that 
are already being experimentally investigated ll27[ [281. We 
propose that experiments be carried out in the SC-graphene 
system to directly confirm our quantum phase transition and 
current switching predictions. 

Acknowledgments. We thank M. Feigel'man, E. Hwang, 
J. Lau and S. Tewari for stimulating discussions. VG acknowl- 
edges the hospitality of Boston University visitors program. 
This work has been supported (RML and SDS) by US-ONR 
and LPS-NSA-CMTC. 



[1 

[2 
[3 
[4 
[5 
[6 
[7 
[8 
[9 

[10 
[11 

[12 
[13 

[14 

[15 
[16 

[17 
[18 

[19 
[20 
[21 

[22; 

[23 
[24; 

[25 
[26 
[27 
[28 



A. M. Goldman and N. Markovic, Physics Today 51, 39 (1998). 
I. S. Beloborodov et al., Rev. Mod. Phys. 79, 469 (2007). 

B. G. Orr et al., Phys. Rev. Lett. 56, 378 (1986). 

H. S. J. van der Zant et al., Phys. Rev. B 54, 10081 (1996). 
Y. Takahide et al., Phys. Rev. Lett. 85, 1974 (2000). 

A. J. Rimberg et al., Phys. Rev. Lett. 78, 2632 (1997). 

G. Schoen and A. D. Zaikin, Phys. Rep. 198, 237 (1990). 
S. Chakravarty et al., Phys. Rev. Lett. 56, 2303 (1986). 

M. V. Feigel'man and A. I. Larkin, Chem. Phys. 235, 107 
(1998). 

V. Galitski and A. I. Larkin, Phys. Rev. Lett. 87, 087001 (2001). 

F. W. J. Hekking and Yu. V. Nazarov, Phys. Rev. B 49, 6847 
(1994). 

C. Bruder, R. Fazio and G. Schoen, Physica B 203 240 (1994). 
R. M. Lutchyn and L. I. Glazman, Phys. Rev. B 76, 104507 
(2007). 

I. L. Aleiner and P. W. Brouwer and L. I. Glazman, Physics 
Reports 358, 309 (2002). 

J. Gonzalez and E. Perfetto, Phys. Rev. B 76, 155404 (2007). 
S. Sachdev, Quantum Phase Transitions (Cambridge University 
Press, 1999). 

J. M. Kosterlitz, Phys. Rev. Lett. 37, 1577 (1976). 

E. Brezin, J. Zinn- Justin and J. C. Le Guillou, J. Phys. A 9, 

L119 (1976). 

W. Hofstetter and W. Zwerger, Eur. Phys. J. B 5, 751 (1998). 

H. Spohn and W. Zwerger, J. Stat. Phys. 94, 1037 (1999). 
S.V. Panyukov, A.D. Zaikin, Phys. Rev. Lett. 67 3168 (1991). 

G. Falci, G. Schoen, and G.T. Zimanyi, Phys. Rev. Lett. 74, 
3257 (1995). 

S. L. Lukyanov and P. Werner, J. Stat. Mech. 11, 11002 (2006). 
J. S. Meyer, A. Kamenev, and L. I. Glazman, Phys. Rev. B 70, 
45310(2004). 

Gil Refael et al., Phys. Rev. B 75, 014522 (2007). 

S. Tewari and J. Toner, Europhys. Lett., 74, 341 (2006). 

H. B. Heersche et al, Nature 446, 56 (2007). 
Xu Du et al., arXiv.org:0710.4984. 



